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Abstract 

We calculate the A-dependence of open charm and beauty as well as Drell- 
Yan pair production cross sections on nuclear targets taking into account the 
difference of quark and gluon distributions in free and bound nucleons. The 
results are presented for both proton- nucleus and nucleus- nucleus collisions. For 
the case of heavy flavour production at comparatively low energies, if a ~ A a , the 
values of a are slightly higher than unity. With increasing energy a decreases 
in all cases and becomes smaller than unity. We also show a as a function 
of Feynman-j; of the produced QQ or l + l~ pair and obtain that it decreases 
significantly in the beam fragmentation region. The sensitivity of our results to 
the used set of parton distributions is also discussed. 
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1 INTRODUCTION 

Processes of heavy quark and lepton pair production on nucleon and nuclear targets at 
high energies are very interesting from both theoretical and practical reasons. These 
processes provide a method to study the internal structure of hadrons in the very small 
x-region. Some limits on the possible non-perturbative contributions to this region can 
also be obtained. Realistic estimates of the cross sections of heavy quark production 
are necessary in order to plan experiments on the existing and future accelerators. 

Predictions are usually obtained in the framework of perturbative QCD in the lead- 
ing and the next-to-leading order a s -expansion. In the case of a nuclear target it is 
usually assumed that the cross section of heavy flavour production (or any other hard 
process), <j(QQ), should be proportional to A a , with a — 1; that is, in agreement with 
the most accurate experimental results of [1,2], where values of a = 1.00 ± 0.05 for all 
.D-meson production by pions at ^Js = 22 GeV and a = 1.02 ± 0.03 ± 0.02 for neutral 
.D-meson production by protons at -\/s = 39 GeV were obtained, respectively. 

However it seems important to consider the problem in more detail because, for ex- 
ample, many experimental results are obtained by extrapolation of results on a nuclear 
target to a nucleon target. It is experimentally well-known [3-9] that parton distribu- 
tions in nuclei (i.e., in bound nucleons) are slightly different from the same distributions 
in free nucleons. So the value of a can differ from unity and it is interesting to estimate 
the size of this effect. Also the differences in the distributions of gluons, sea quarks and 
valence quarks are not the same [|1(J. Moreover, in different hard processes different 
QCD subprocesses give the main contributions, so the effects on the a-values can also 
be different. 

In this paper we calculate the values of a as a function of energy and Feynman-x 
for the cases of open charm, beauty and different mass Drell-Yan pair production in 
proton-nucleus and nucleus-nucleus collisions. 

In our previous paper |TI[] we showed that the results for the A-dependence of heavy 
flavour production depend very weakly on the used set of parton distributions. Three 
different sets, namely MRS-1 JIJ, MT S-DIS [lj and GRV HO [fLj, which can be 



found in CERN PDFLIB [IS|, gave practically the same results for the a behaviour 



(see Table 1 of Ref. [|ll|]). Here we present the results for two sets, GRV HO Jl4 
and MRS-A [Hjj which are in comparatively good agreement with the new HERA data 
and we see that they give practically identical results for both open heavy flavour and 
Drell-Yan pair production in proton-nucleus and nucleus- nucleus collisions. 
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2 CROSS SECTIONS OF HARD PROCESSES IN 
QCD 



The standard QCD expression for the cross section of heavy quark production in a 
hadron 1 - hadron 2 collision has the form 
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a ab - QQ (s,m Q ,^), (1) 



where x a o = 4mg/s and x b o = ^thq/ (sx a ). Here G a /x(x a , /i 2 ) and G b /^(x b , /i 2 ) are 
the structure functions of partons a and b inside hadrons 1 and 2 respectively, and 
s, rriQ, fi 2 ) is the cross section of the subprocess ab — > QQ as given by standard 
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ab^QQ , 



QCD. The latter depends on the parton center-of-mass energy s = (p a + p b ) 2 = x a x b s, 
the mass of the produced heavy quark ttiq and the QCD scale /i 2 . Eq. (]!]) should 
account for all possible subprocesses ab — > QQ. The parton-parton cross section e ab ^QQ 
can be written in the form |L7| 



a 2 s (^ 



-fab(p, fJ 2 , m 2 c 



with 



P = 4m 2 /s 



(2) 



(3) 



and 



fab(p,p 2 ,m 2 Q ) = f^(p)+^a s ^ 2 )[fH\p) + tii\p)\n(v 2 /m 2 Q )} . (4) 

The functions fa b \p), fi b \p) and fi b \p) can be found in |L7j for heavy flavour 
production. 

Formulae with the same structure as (l)-(4) can be found in |L8| for Drell-Yan pair 
production. In this case M t +i- plays the role of 2m<g in x a o, x b Q and cr ab ~^ QQ . 

For the numerical calculations we wrote the nuclear structure function G b /^(x b , p 2 ) 
in the form 



G b/A (x b , jj 2 ) = A - G b/N (x b , jj 2 ) ■ R b (x, /i 2 ), (5) 

similarly to Ref. |0| and we take the values of i?^(x,/i 2 ) for gluons, valence and sea 
quarks from Ref. [10]. They are presented in Fig. 1. The values of R^(x, fx 2 ) in [10] are 
given for x > 10 -3 which is not small enough at high energies. So at x < 10~ 3 we used 
two variants of the R^(x, fi 2 ) behaviour for gluon and sea quark distributions: The first 
is the constant frozen at x — 10~ 3 (dashed-dotted curves for gluons and solid curves 
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for sea quarks in Fig. 1). The second is the extrapolation as x^ (dotted curves in Fig. 
1) of the distribution which gives the main contribution to our cross section, with (3 = 
0.096 and 0.040 for charm and beauty production respectively (only for gluons), and 
(3 = 0.109, 0.096 and 0.072 for Drell-Yan pair invariant masses M 2 +l _ = 5, 25 and 100 
GeV 2 respectively (only for sea quarks). Such behaviour is in qualitative agreement 
with the results of [ f2l| . 

In the case of nucleus-nucleus collisions the parton distributions in both incident 
nuclei should be written in the form (5). It is necessary to note that in such a way the 
main part of the shadowing processes is taken into account. 

3 A-DEPENDENCE OF HEAVY FLAVOUR PRO- 
DUCTION 

In the case of charm production we have used the values m c = 1.5 GeV and fi 2 = 4 
GeV 2 {jj 2 = 5 GeV 2 for the MRS-A set) and in the case of beauty production mj = 5 
GeV and fi 2 = m 2 . 

The obtained results for a determined from the ratios of heavy quark production 
cross section on a gold target and on the proton are presented in Fig. 2 for the GRV 
HO and MRS-A sets[] and the two variants of gluon distribution ratios in the small 
x-region. Here ^s^n is the cm. energy for the interaction of the incident proton 
with one target nucleon. One can see that at fixed target energies the values of a 
are slightly higher than unity, which is not in contradiction with the results of Refs. 
[1,2]. However a decreases with increasing energy and this effect is larger in the case 
of charm production than in the case of beauty. One can see also that the difference 
between the two variants for i?^(x,/i 2 ) at x < 10~ 3 becomes important only at the 
highest energies. 

We also calculate the values of a for different Feynman-rc (xf) regions using xf = 
x a — x b at energies ^snn = 39 GeV and 1800 GeV. The results are presented in Fig. 
3. At negative and moderate xf (in the nucleus fragmentation region) the values of a 
are slightly higher than unity. However in the case of charm production in the beam 
fragmentation region (positive xp) the values of a become essentially smaller than 
unity. For beauty production the last effect is expected only at very high energies. 

The obtained results for A-dependences of charm and beauty production in the 
symmetric case of gold-gold collisions are presented in Figs. 4 (a as a function of 
initial energy) and 5 (a as a function of xp at two energies). Here ^/snn is the cm. 

all the calculations in this paper little difference, if any, has been found in the results for a 
between these two sets of parton distributions. 
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energy for one nucleon-nucleon interaction. One can see a qualitative agreement with 
the case of proton-nucleus collisions (with the difference that the trivial value of a is 
here equal to two). Again we can see that the values of a are dependent on the initial 
energy and Xp if the energy is high enough. 

Calculations for charm production in perturbative QCD (including preequilibrium 
charm production from secondary minijet gluons) can be found in |21| for central gold 
on gold collisions. 



4 A-DEPENDENCE OF DRELL-YAN PAIR PRO- 
DUCTION 

As said above, the QCD expression for heavy lepton pair production in a hadron 1 
- hadron 2 collision has the same form as Eq. (1) but with another matrix element 
which can be found in Ref . |L8) . Now we have an additional variable - the mass of the 



produced pair M\+\- which plays more or less the same role as the mass of the heavy 
quark. However now its value can be measured experimentally and do not lead to any 
uncertainty. We have used the values of QCD scale fi 2 = Mf +l -. 

The obtained results for a determined in the same way as in heavy quark production, 
i.e., from the ratios of Drell-Yan pair production cross section in proton-gold and 
proton-proton collisions, are presented in Fig. 6 for the GRV HO and MRS-A sets and 
the two variants of sea quark distribution ratios in the small x-regionf]. Here again 
y/ s NN is the cm. energy for the interaction of the incident proton with one target 
nucleon. Contrary to the case of heavy flavour production the values of a are never 
higher than unity. They decrease more or less monotonically with increasing energy 
and this effect becomes smaller with increasing mass. Again the two sets of parton 
distributions give practically identical results for the energy dependence of a and the 
difference between the two variants for R^(x,fi 2 ) at x < 10~ 3 becomes important 
only at the highest energies. The predicted Mi+i- dependence of a for Drell-Yan pair 
production in proton-nucleus collisions at different energies is also shown in Fig. 6. 

The results of our calculations of a as a function of xp for heavy lepton pair pro- 
duction are presented in Fig. 7. The qualitative picture is similar to the case of heavy 
flavour production; the numerical difference of our predictions from the value a = 1 is 
here more significant for the case of a not very large mass of the lepton pair. 

Predictions of A-dependence for Drell-Yan pair production cross section in symmet- 
ric nucleus- nucleus collisions are shown in Figs. 8 and 9. As for the case of proton- 



2 As is well-known gluons are dominant for heavy flavour production, while sea quarks dominate 
Drell-Yan production. 
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nucleus collisions, all the effects are qualitatively the same but numerically larger than 
in the case of heavy flavour production. Our results are in agreement with the calcu- 
lations of 151. 



5 CONCLUSIONS 

In this paper we calculate the A-dependence of charm and beauty as well as Drell- 
Yan pair production using standard QCD formulas and accounting for the difference 
of parton distributions in free and bound nucleons. If one parametrize these cross 
sections as a ~ A a , the value of a is slightly different from unity at the available 
energies. For the case of heavy flavour production at comparatively low energies the 
obtained values of a are a little higher than unity. This should be connected with some 
nucleon-nucleon correlations which change the large-x parton distributions. In the case 
of Drell-Yan pair production such effect is absent because of the different relations in 
the contributions of valence quarks, sea quarks and gluons. 

At higher energies the values of a decrease and become smaller than unity. At 
y/SNN = 1800 GeV we expect a value of a ~ 0.95 for charm and low-mass Drell-Yan 
pair production. The decrease of the ratio R^(x,n 2 ), which results in a decrease of 



a, can be connected with the effects of parton density saturation |23| which in heavy 
nuclei occur at x-values higher than in the proton. 

If we consider two small and different values of x a and Xf, in Eq. (1) for the case of 
proton-nucleus interaction, it is clear that the contribution to the inclusive cross section 
from the region x a < Xb should be larger than the mirror contribution (x a > Xj,) because 
the value of the ratio Rj^(x, /i 2 ) in the first case is larger. It means that heavy quark and 
Drell-Yan pairs will be produced preferably in the nucleus fragmentation hemisphere, 
i.e., asymmetrically, which is quite usual and has been confirmed experimentally in the 
case of light quark production. From Figs. 3 and 7 it is clear that charm and low- 
mass Drell-Yan pair production on nuclear targets at LHC energies will give important 
information on the nuclear shadowing of the structure functions at small x. 

Besides, the experimental measurements of the effects predicted for different hard 
interactions should allow (in principle) to control the validity of the conventional ex- 
traction of parton distributions from the experimental DIS data. 

Let us repeat again that almost all nuclear shadowing corrections are accounted for 
in our calculations because they contribute to the Rb(x,fi 2 ) ratios, which have been 
extracted from the experimental data. On the other hand these shadowing corrections 
are not numerically large in the considered hard processes. So we can assume that 
the corrections which are not accounted for do not change significantly the obtained 
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results. 
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Figure captions 



Fig. 1. Functions i?^(a;,/i 2 ), Ry(x,fi 2 ) and Rg(x,[i 2 ), which determine the ratios 
of the distributions for protons in the nucleus versus free protons, for gluons (dashed- 
dotted and dotted curves, see text), valence quarks (dashed curves) and sea quarks 
(solid and dotted curves, see text) respectively, for /i 2 = 5 GeV 2 (upper figure) and 
fj, 2 = ml (lower figure). 

Fig. 2. Energy dependence of a for charm and beauty production in proton-gold 
collisions for GRV HO (solid and dashed curves) and MRS-A (dotted and dashed- 
dotted curves) structure functions and using extrapolated (dashed and dashed-dotted 
curves) and frozen at x — 1CT 3 (solid and dotted curves) ratios of gluon distributions. 
Fig. 3. Feynman-x dependence of a for charm and beauty production in proton-gold 
collisions at v /sjvjv = 39 GeV (upper figure) and 1800 GeV (lower figure) for GRV HO 
and MSR-A structure functions and using extrapolated and frozen at x = 10~ 3 ratios 
of gluon distributions (with the same conventions as in Fig. 2). 

Fig. 4. Energy dependence of a for charm and beauty production in gold-gold colli- 
sions for GRV HO and MRS-A structure functions and using extrapolated and frozen 
at x — 10~ 3 ratios of gluon distributions (with the same conventions as in Fig. 2). 
Fig. 5. Feynman-x dependence of a values for charm and beauty production in gold- 
gold collisions at y/s NN = 39 GeV (upper figure) and 1800 GeV (lower figure) for GRV 
HO and MSR-A structure functions and using extrapolated and frozen at x — 10~ 3 
ratios of gluon distributions (with the same conventions as in Fig. 2). 
Fig. 6. Mass (upper figure) and energy (lower figure) dependence of a for Drell-Yan 
pair production in proton-gold collisions for GRV HO and MRS-A structure functions 
and using extrapolated and frozen at x — 10~ 3 ratios of sea quark distributions (with 
the same conventions as in Fig. 2). 

Fig. 7. Feynman-a; dependence of a values for heavy lepton pair production in proton- 
gold collisions at ^/snn = 39 GeV (upper curves in each figure) and 1800 GeV (lower 
curves in each figure) and different masses for GRV HO and MSR-A structure functions 
and using extrapolated and frozen at x — 10~ 3 ratios of sea quark distributions (with 
the same conventions as in Fig. 2). Note that at y / sjv r /v — 39 GeV all curves coincide. 
Fig. 8. Mass (upper figure) and energy (lower figure) dependence of a for Drell-Yan 
pair production in gold-gold collisions for GRV HO and MRS-A structure functions 
and using extrapolated and frozen at x — 10 -3 ratios of sea quark distributions (with 
the same conventions as in Fig. 2). 

Fig. 9. Feynman-rr dependence of a values for heavy lepton pair production in gold- 
gold collisions at ^s^n — 39 GeV (upper curves in each figure) and 1800 GeV (lower 
curves in each figure) and different masses for GRV HO and MSR-A structure functions 
and using extrapolated and frozen at x — 10~ 3 ratios of sea quark distributions (with 



8 



the same conventions as in Fig. 2). Note that at ^/snn = 39 GeV all curves coincide. 
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